Critical generalized inverse participation ratio distributions 
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The system size dependence of the fluctuations in generalized inverse participation ratios (IPR's) 
I a (q) at criticality is investigated numerically. We focus on a three-dimensional (3D) system with 
unitary symmetry, a 2D system with symplectic symmetry and a ID system with orthogonal sym- 
metry. The variances of the IPR logarithms are found to be scale-invariant at the macroscopic 
limit. The finite size corrections to the variances decay algebraically with nontrivial exponents, 
which depend on the Hamiltonian symmetry and the dimensionality. The large-q dependence of the 
asymptotic values of the variances behaves as q 2 according to theoretical estimates. These results 
ensure the self-averaging of the corresponding generalized dimensions. 
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The statistical properties of wave functions have been 
intensively studied during recent years, helping in our un- 



lena in different areas of physics, 
atomico to microwaves and 




derstanding of phen< 
ranging from nucle 

mesoscopic physics. Of considerable interest is 
what happens at the critical point of a metal-insulator 
transition (MIT), in which eigenstates presenting strong 
amplitude fluctuations were demonstrated to be multi- 
fractal objects. 

The Anderson MIT essentially depends on the dimen- 
sionality and symmpi|ries of the system. The field the- 
oretical formulationliy of the MIT shows that the crit- 
ical behavior can be described within a framework of 
three universality classes: orthogonal (/3 = 1), unitary 
(P = 2), and symplectic [fi — 4). The symmetry pa- 
rameter (3 is the number of independent real components 
that characterizes a matrix element of the corresponding 
Hamiltonian. The two major symmetries in field theory 
are time-reversal symmetry and spin-rotation symmetry. 
The system belongs to the orthogonal class if it has both 
symmetries, to the unitary class if time-reversal symme- 
try is broken, and to the symplectic class if the system has 
time-reversal symmetry but spin-rotation is broken. The 
relevant terms in the hamiltonian are a coupling to an 
applied magnetic field, which breaks time-reversal sym- 
metry, and the spin-orbit interaction, which breaks spin- 
rotation symmetry. 

Castellani and PelitiEiJ considered the d-^= 2 + e (e <C 1) 
expansion of Wegncr's nonlinear a modeE3 for the gener- 
alized inverse participation ratios (IPR's) of eigenstates, 



I a (q) = f d d r \Mr)\ 2q 



J^-d q (q- 



(1) 



and concluded that the wave functions ip a (r) at the tran- 
sition point of a MIT show multifractal behavior charac- 
terized by an infinite set of generalized fractal dimensions 
d q . The index a labels the different cigenfunctions and f2 
denotes a d-dimensional region with linear dimension L. 
Equation ([!]) is valid for individual states and for their 
ensemble average since the spectrum of multifractal di- 
mensions has universal features for states in the vicinity 



of the MIT0. 

As regards the IPR fluctuations, we shall first review 
the existing analytical results. When these fluctuations 
were studied for two-dimensional systems, in the frame- 
work of the supersymmetry methodJjoEj it was found 
that the distribution function of I a (q), normalized to its 
typical value, is scale-invariant at criticality.O Although 
this case does not present a true Anderson transition, 
the above result motivated the conjecture that in general 
the normalized distribution function of I a (q) is univer- 
sal, i.e., size independent for L — > oo. Accordingly, it is 
assumed that the distribution function of In I a (q) is a uni- 
versal curve that, while horizontally shifted by changes 
in L, keeps the same, form. Within the same approach, 
Mirlin and EversE3E3 study the fractal properties of the 
power-law random banded matrix (PRBM) ensemble at 
criticality. This model is characterized by a parameter, b 
[see Eq. (^) below]. In the two limiting cases, b^$> I and 
b <c 1, for which analytical treatment is feasible, they 
found the distribution function of the IPR to be scale- 
independent at the macroscopic limit. Particularly, they 
analyzed the behavior of the variance of the distribu- 
tion function of I a (q), normalized to its typical value, in 
the limit b <C 1 . For values of q larger than some critical 
value, q c « 2.41, where the IPR distribution is dominated 
by its slowly decaying power law tail, the same authors 
obtained for L — > oo 



(oo) = (q/q c ) 



(2) 



What we are trying to point out is that a general ana- 
lytical description of the statistical properties of the gen- 
eralized IPR at the critical point of realistic MIT's (con- 
ventional Anderson transition, quantum Hall transition, 
transition in d = 2 for electrons with strong spin-orbit 
coupling, etc.) is still lacking. This is why we addressed 
the problem using numerical calculations. To our knowl- 
edge, no such calculations have been previously reported. 

Our aim was to find the system size dependence of the 
generalized IPR distributions in order to check whether 
they converge at large L to the conjectured scale invari- 
ant distributions or not. A second motivation for this 
work was to clarify whether the corresponding general- 
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ized dimensions are well-defined in the macroscopic limit 
or presents scale-invariant distributions. 

Using the exact eigenstates from numerical diagonal- 
izations, we obtained the size dependence of the distri- 
bution function Vp(z q = hil a (q)) at the critical point 
of a three-dimensional (3D) system with unitary sym- 
metry, a two-dimensional (2D) system with symplectic 
symmetry and a one-dimensional (ID) system with or- 
thogonal symmetry. All these systems possess a mobility 
edge and are modeled by the corresponding tight-binding 
Hamiltonian. We found that as L increases Vp becomes 
wider and the height of its peak smaller with a tendency 
towards saturation. Therefore, their variance increases 
as L increases and will eventually saturate to a constant 
finite value for L — > oo. In order to extrapolate to macro- 
scopic systems, we propose a finite size correction to the 
variance a 2 of z q according to 
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with a~ (oo), a q and y q > being three adjustable pa- 
rameters. Practically for all triads {q, (3, d} we found that 
the exponent y q was always close to the generalized di- 
mension d q divided by 2j3d. Thus, for the exponent y q 
we propose the following equation, 



FIG. 1: Finite size correction S 2 z q (L) as a function of system 
size L for the 3D disordered system with unitary symmetry 
and several q values: q — 2 (circles), q — 4 (squares), q = 6 
(diamonds), and q = 8 (triangles). The straight lines are 
linear fits to Eq. (ta). 
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u 1 
2/3d 



(4) 



and keep only two free parameters. This relation con- 
stitutes the main result of the present publication. As 
we will see below, our numerical data strongly support 
this conjecture. Equation (0) is a generalisation of the 
exponent 7 = d%/2d found by our groupcJ for I a {2) in 
disordered systems with orthogonal symmetry. 

Some technical details follow: the critical eigenfunc- 
tions (and eigenvalues) of the Hamiltonian matrices, Eqs. 
@> (§)' an d ©1 are obtained by numerical diagonal- 
ization. In 3D and 2D systems we use techniques for 
large sparse matrices Jla while for the ID case standard 
diagonalization subroutines are used, since we must deal 
with full matrices. We consider a small energy window 
(— T), 77) around the center of the spectral band, where 
r] = 1 (ry = 0.4) for 3D and 2D (ID) systems. Reduc- 
ing the width of the above windows does not alter the 
results. For an accurate extrapolation to the thermo- 
dynamic limit, the system size changes by one order of 
magnitude for 3D and 2D systems and by two orders in 
the case of ID. The number of random realizations is such 
that the number of critical states included for each L is 
roughly 1.5 x 10 5 , while, in order to reduce edge effects, 
we impose periodic boundary conditions in all the cases 
considered. WeJiave checked that the normalized nearest 
level variancesEj are indeed scale-invariant at each criti- 
cal point studied. Below we present our numerical results 
and compare them to the conjecture (f|). 

3D unitary symmetry. First we focus on the breaking 
of time-reversal symmetry by a constant applied mag- 
netic field. We consider the standard Anderson model 



with diagonal disorder and Peierls phase factors in the 
hopping matrix elements 



(5) 



where c' n (c„) are the creation (annihilation) operators of 
an electron at a site n of a simple-cubic lattice, and the 
summation in the second term is over nearest neighbors. 
The diagonal energies are randomly distributed with con- 
stant probability in the interval —W/2 < e n < W/2, 
where W parametrizes the strength of random poten- 
tial. Energies are measured in units of the modulus 
of the hopping matrix elements, V nm , and lengths in 
units of the lattice constant. The hopping matrix ele- 
ments, V nm = e T2 " az for neighbors in the y direction 
and V nm — 1 in any other direction, describe a system 
with a homogeneous magnetic field, B, in the x direc- 
tion, while the Peierls phase a = eB/hc is the number of 
flux quanta <&o = hc/e threading a lattice cell. Here, the 
Landau gauge with the vector potential A = (0, — Bz, 0) 
is chosen. The system size varies between 4 and 30, and 
for the critical disorder we take the value W c — 18.35 
reported in Ref. 20 for a = 0.25. 

In Fig. 1 we represent on a log-log scale the finite size 
correction 5 2 z q (L) = a 1 (00) — a 1 (L) as a function of 
L for the 3D disordered system with unitary symmetry 
and different values of q: 2 (circles), 4 (squares), 6 (di- 
amonds), and 8 (triangles). The good agreement with 
Eqs. (g) and (|J) is evident. The fitted values of the 
free parameters are a-i = 2.01, 04 = 19.71, a$ — 52.88, 
as = 100.7, while the corresponding a z (00) are summa- 
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rized in Fig. 4 (squares). The slopes of the straight lines 
have not been fitted and correspond to Eq. (J4|), with 
= 2, d = 3, d 2 = 1.37, di = 1.07, d 6 = 0.98, and 
c?8 = 0.94. These generalized dimensions were estimated 
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from the slope, d q (l 



■ q), of {z q ) versus InL. In order to 
produce a diagram of reasonable size the curves for q = 2, 
4, and 6 have been shifted vertically. The small difference 
in the slopes for q > 4 correspond to the slight differences 
in the values of d q . The value a Z2 (oo) = 1.40 found is in 
good agreement with the conjecture a Z2 (oo) ~ 1.0 

2D symplectic symmetry. In order to calculate the crit- 
ical eigenstates (two-component spinors) Ando's model 
HamiltonianEil was used, 



where c], a (c n ,a) are the creation (annihilation) opera- 
tors of an electron at a site n of a square lattice with 
the spin component a and m denotes the sites adjacent 
to site n. The on-site energy e n is randomly distributed 
around zero according to a box distribution with a width 
W, which specifies the degree of disorder. The index a 
takes on values of +1 and —1, and denotes spin- up or 
spin-down. The transfer matrices V n>CT - >m<CT i = V x or V y 
depend on the direction of the nearest neighbor site 



V x = 



Vx V 2 
-V 2 Vx 



Vy = 



Vi ~iV 2 
-iV 2 Vx 



and the strength of the spin-orbit coupling is given by the 
parameter S = V 2 /V, with V = (Vf +F 2 2 ) 1/2 taken to be 
the unit of energy. We chose S =1/2 and, for the critical 
disorder, W c = 5.98, which was previously found in Ref. 
22. The system size ranges from L = 10 to L = 120. 

Figure 2 shows the same size corrections as Fig. 1 
for the 2D disordered system with symplectic symmetry. 
In this case the fitted values of the free parameters are 
a 2 = 0.17, 04 = 5.10, clq — 15.77, and a§ = 31.27 and the 
values of a 2 z (oo) are shown in Fig. 4 (circles). As in the 
3D case, the slopes of the straight lines correspond to Eq. 
(@), with (3 = 4, d = 2, d 2 = 1.60, d A = 1.34, d 6 = 1.23, 
and d$ = 1.18. These dimensions were also estimated in 
a similar way. Again, the value found, a Z2 (oo\ = 0.46, is 
in agreement with the conjecture a Z2 (oo) ~ 1.0 

ID nrihnrjnunl symmetry. In this case the PRBM 
modelli3ii3 is used. The Hamiltonian, which describes a 
disordered ID sample with random long-range hopping, 
is represented by real symmetric matrices, whose entries 
are randomly drawn from a normal distribution with zero 
mean (7i nm ) = 0, and a variance depending on the dis- 
tance between lattice sites 



(\H nm \ 2 ) - 1 + ( | n _ m | /6 j 2 



(7) 



The model describes a whole family of critical theories 
parametrized by b, which determines the critical dimen- 
sionless conductance, in the same way as the dimension- 
ality labels the different Anderson transitions. In the 
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FIG. 2: As for Fig. 1, for the 2D disordered system with 
symplectic symmetry. 
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FIG. 3: As for Fig. 1, 
orthogonal symmetry. 



for the ID disordered system with 



intermediate regime b ~ 1, where there is no analytical 
solution, we obtained two different diverging exponents 
for the [Correlation and localization lengths at the critical 
region. E3 The system size ranges between L = 100 and 
L = 10 000, and for the calculations we choose b = 0.2. 

Figure 3 shows the same size corrections as Fig. 1 
for the ID disordered system with orthogonal symmetry. 
The fitted values of the free parameters are a 2 = 0.47, 
04 = 4.28, ag = 11.16, a§ = 20.95, and cr 2 (oo) are shown 
in Fig. 4 (diamonds). As in the two previous cases, the 
slopes of the straight lines correspond to Eq. (0), with 
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FIG. 4: Asymptotic values of the standard deviation a Zq (oo) 
as a function of q for the ID system with orthogonal sym- 
metry (diamonds), 2D symplectic (circles), and 3D unitary 
(squares). 



(3 = 1, d = 1, d 2 = 0.36, d A = 0.29, d 6 = 0.27, and 
d$ = 0.26. The last four dimensions were estimated in 
a way similar to those for Fig. 1. As in the 3D and 2D 
cases the value a Z2 (oo) = 0ij7 found is in agreement with 
the conjecture a Z2 (oo) ~ 1.0 

We checked that Eqs. (^) and (^) adequately describe 
the IPR fluctuations in the standard 3D Anderson tran- 
sition (orthogonal symmetry) and in the critical point 
of the ID system with long-range hopping, Eq. (Q), 
when unitary and symplectic symmetries are included. 
We are confident that further analytical development in 
IPR statistics at criticality will confirm our conjectured 
result. 

The statistical properties of wave functions are inti- 
mately related to the universal conductance fluctuations. 
As was shown in Ref. 6 , the relative variance of I a (2) is 
of the order of 1/g 2 , g being the dimensionless conduc- 
tance (in units of e 2 /h). The extrapolation of this rela- 
tion to the transition point, in which g ~ 1, is the most 
direct connection between our results and experiments. 
In this way, our results could be checked in experiments 
with random systems (semiconductor heterostructures, 
metal oxide semiconductor MOS inversion layers, vapor- 
deposited films, doped and amorphous semiconductor, 



etc.) using samples with varying size. 

Finally, we focus on the q dependence of the asymp- 
totic values a Zq (oo) for each universality class. In Fig. 4 
we represent the standard deviation o~ Zq (oo) versus q for 
the three transitions studied: ID orthogonal (diamonds), 
2D symplectic (circles), and 3D unitary (squares). The 
slopes of the fitted straight lines are 0.63, 0.86, and 1.39, 
respectively. It is interesting to note that the nonlinear 
a model estimates, Eq. (0), of the ID PRBM model is 
valid for all three transitions. These results further sup- 
port our conjecture, described by Eq. (Q). The different 
values of the slopes directly reflect the universality class 
and the dimensionality of the system. From the above 
slope for the ID system we obtain q c = 1.59, which is 
less than the value of 2.41 predicted theoretically. The 
origin of the deviation is that the latter value of q c was 
obtained for the distribution function of I a (q) normal- 
ized to its typical value, while the former q c corresponds 
to the distribution function of In I a (q) . Anyway, we have 
estimated the value of q c from the extrapolated long tails 
of the normalized I a {q) distributions. The obtained value 
q c = 2.21 is in close agreement with the predicted one for 
these distributions. 

If we now assume that the fluctuations of the gener- 
alized dimensions d q are due-to fluctuations of the cor- 
responding generalized IPR,o from Eq. ([!]) we obtain 
a\(V) = a Zq (L)/(q — l) 2 ln 2 L for the variance of d q . 
From this relation, and taking into account Eqs. (Q), ||), 
and (^) , we see that a Zq (L) vanishes at a rate 1/ In L when 
L — > oo. Therefore, we conclude that all generalized di- 
mensions are self-averaged quantities, i.e., they tend to 
well-defined single values at the macroscopic limit. 

In summary, we have performed a detailed numerical 
analysis of the statistical properties of the IPR distribu- 
tions at the critical point of three standard MIT's, find- 
ing the finite asymptotic values of the variances of the 
distributions in each case. The corresponding finite size 
corrections decay with system size as a power law with 
exponents that depend on the system dimensionality and 
on the Hamiltonian symmetry. For large values of q, the 
asymptotic values a Zq (oo) oc q 2 for all three transitions 
is in agreement with the analytical estimates of the ID 
PRBM model. Therefore, all generalized dimensions are 
well-defined at the macroscopic limit. 
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